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Abstract 

An analytical model of the electrostatic force between the tip of a non-contact Atomic Force 
Microscope (nc-AFM) and the (001) surface of an ionic crystal is reported. The model is able 
to account for the atomic contrast of the local contact potential difference (CPD) observed while 
nc-AFM-based Kelvin Probe Force Microscopy (KPFM) experiments. With the goal in mind to 
put in evidence this short-range electrostatic force, the Madelung potential arising at the surface 
of the ionic crystal is primarily derived. The expression of the force which is deduced can be split 
into two major contributions: the first stands for the coupling between the microscopic structure of 
the tip apex and the capacitor formed between the tip, the ionic crystal and the counter-electrode; 
the second term depicts the influence of the Madelung surface potential on the mesoscopic part of 
the tip, independently from its microscopic structure. The former has the lateral periodicity of the 
Madelung surface potential whereas the latter only acts as a static component, which shifts the 
total force. These short-range electrostatic forces are in the range of ten pico-Newtons. Beyond 
the dielectric properties of the crystal, a major effect which is responsible for the atomic contrast 
of the KPFM signal is the ionic polarization of the sample due to the influence of the tip/counter- 
electrode capacitor. When explicitly considering the crystal polarization, an analytical expression 
of the bias voltage to be applied on the tip to compensate for the local CPD, i.e. to cancel the 
short-range electrostatic force, is derived. The compensated CPD has the lateral periodicity of 
the Madelung surface potential. However, the strong dependence on the tip geometry, the applied 
modulation voltage as well as the tip-sample distance, which can even lead to an overestimation 
of the real surface potential, makes quantitative KPFM measurements of the local CPD extremely 
difficult. 

PACS numbers: 07.79.Lh, 41.20.Cv, 73.40.Cg 
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I. INTRODUCTION 

Electrostatic forces play a key role in non-contact Atomic Force Microscopy (nc-AFM), 
not only in the imaging process^ but also for the investigation of the surface electronic prop- 
erties. Electronic properties such as the work function and surface charges can be acquired 
by Kelvin Probe Force Microscopy (KPFM)^'^ simultaneously to imaging topography by 
nc-AFM. In KPFM, a feedback is used to apply a voltage between the tip and the sample in 
order to minimize the electrostatic tip-sample interaction. For metals, this voltage is equal 
to the contact potential difference (CPD), i.e. the bias voltage to be applied between the tip 
and the surface to align their fermi levels. It is connected to the difference between the work 
functions of the two surfaces, and thereby to their local electronic properties, according to: 

^0 = (Ptip — (Psample = QVcpd, (l) 

q being the elementary electrical charge: q = 1.6 x 10^^^ C 

Nowadays, two KPFM-based techniques provide facilities to map the spatial variations of 
the CPD on the nanometer scale, namely Frequency-^ or Amplitude- Modulation-KPFM^i^ 
(FM- or AM-KPFM, respectively). These methods were primarily applied to metallic and 
semiconducting surfaces to study the distribution of dopants in semiconductors'^, or the 
adsorption of organic molecules (for an overview see ref.|9|]). In a few experiments, even 
molecular— or atomioiiii^ii^ contrast has been reported. The extension of the technique to 
insulating surfaces, was performed more recently, as demonstrated by the results reported 
on thin ionic films on metals^*^, or on the contribution of bulk defects to the surface charge 
state of ionic crystals^*^. 

In this work, atomic corrugation of the CPD signal is reported for the first time on the 
(001) surface of a bulk ionic crystal of KBr. For that purpose, KPFM experiments were 
carried out in ultrahigh vacuum with a base pressure below 10^^° mbar using a home built 
non-contact atomic force microscope operated at room temperature^. A highly doped silicon 
cantilever with a resonance frequency /o ~ 160 kHz, and a spring constant k ^ 21 N.m~^ 
was used. The typical oscillation amplitude of the fundamental bending resonance was 
^ 10 nm. The cantilever was annealed (30 min @ 120°) and gently sputtered with Ar"*" ions 
(1-2 min @ 680 eV). The KBr crystal was cleaved in ultrahigh vacuum along the (OOl)-plane 
and subsequently annealed at 120° during an hour. The KPFM signal was detected using the 



AM mode, as described in detail in Ref. |5|. While these experiments, atomic-scale contrast 
was as well visible in the topography channel (data not shown). The CPD measurements are 
reported in figsH^ and b. In figHt, the image exhibits atomic features, the measured period 
of which is 0.63 nm, which is visible in the joint cross section. This value matches to a good 
agreement the lattice constant of KBr, 0.66 nm. The vertical contrast yields about 100 mV 
around an average value of -3.9 V, the origin of which will be discussed in section HVBI A 
striking aspect of those results is the robustness of the imaging process in terms of stability 
and reproducibility upon the tips used. These results suggest an intrinsic imaging process 
relying on the microscopic origin of the contact potential arising at the sample surface. In 
this case, the CPD rather turns into local CPD, consistently with the concept of local work 
function which has been introduced by Wandelt on metals^^. By "intrinsic imaging process", 
it is meant that the contrast can be accounted for with a tip consisting of a single material, 
namely a metal. Thus, the atomic contrast neither relies on adsorbed nor on unstable species 
at the tip apex, as often reported for topographic^^ or dissipation^i imaging by nc-AFM. 

In this work, in order to understand the local CPD contrast formation, an electrostatic 
model is proposed that allows us to derive an expression of the short-range electrostatic force 
occurring between a biased metallic tip of a nc-AFM microscope and the surface of a bulk 
dielectric. On the contrary to more refined numerical methods proposed in the literature for 
almost similar systems^^i^^, the analytical development is restricted to a simple tip geometry 
and a classical, continuous electrostatic approach. Notwithstanding, the model allows us to 
define a general frame, able to put in relation the surface electrostatic properties with the 
imaging process yielding atomic contrast of the CPD on ionic surfaces. Obviously, the main 
results presented here remain qualitatively correct for more complex tip geometries, although 
numerical methods are then required to get quantitative numbers. 

The motivations for that work are twofold. On the one hand, to our knowledge, a com- 
pact modelization of the short-range electrostatic forces responsible for the atomic contrast 
in KPFM on a bulk ionic crystal is still lacking. On the other hand, when evaluating KPFM 
experiments, the relative complexity of FM- or AM-KPFM experimental setups, both in- 
cluding four electronic controllers, makes the interpretation of the experimental images, and 
primarily CPD images, tedious, especially when dealing with atomic resolution. Several 
groups analyzed the KPFM imaging mechanism in order to evaluate their data in terms of 
quantitative values and lateral resolution^'^'^'^i'^'^^'^'^i^'^, a few of them also compared 



AM- and FM-KPFM in terms of evaluating the force and its gradient, respectively^*^. The 
fact that many of the above mentioned calculations point out that KPFM results, especially 
for nano-objects, show a strong distance dependence, clearly points out that for a careful 
analysis, it is not sufficient to only calculate the electrostatic tip-sample interaction. It is 
rather important to perform simulations including all imaging mechanisms and in particular 
also the distance control in order to exclude artifacts due to the feedback circuits. In order to 
explore the origin of the CPD atomic contrast, we are aiming to closely mimic a real KPFM 
setup by means of an earlier developed nc-AFM simulator—. In the present case of an ionic 
surface, it is required to consider a large slab of ions, virtually infinite, to describe properly 
the electrostatic interaction. This is hardly feasible by means of ab initio calculations which 
fail to describe tip-surface systems involving a too big number of atoms. Therefore, prior to 
simulating the CPD contrast on ionic surfaces by means of our simulator, which is the scope 
of a future work, we have to find an analytical expression for the electrostatic tip-sample 
interaction. 

The following section details the boundary-value electrostatic problem leading to the 
expression of the force (section IIII|) . In section IIVI the analytical expression of the local 
CPD is derived and discussed to be put in relation with the experimental observations. The 
implications for KPFM experiments are discussed as well. 

II. ELECTROSTATIC MODEL 

The geometry of the problem of FM- or AM-KPFM experiments applied to bulk insulating 
materials is reported in fig|2]D. The dielectric sample is an alkali halide crystal like NaCl, 
KBr, KCl with a fee structure and a lattice constant a {cf. fig 12^). In the area where 
the tip is, the model assumes that the crystal carries neither net charge, nor local dipole. 
Its thickness h^ is much larger than all other distances of the problem, typically a few 
millimeters. Below the surface, the crystal is treated as a continuous dielectric medium with 
a dielectric permittivity ea- At the surface, the atomic corrugation of the crystal is described 
by a single layer of alternate point charges arranged with a fee structure perpendicular to 
the [001] direction. The layer extends infinitely in the plane direction. The motivation for 
such a rationalization of the problem will be justified in section ( 1IIB[) . 

The crystal lies on a metallic sample holder (hereafter referred to as the counter- 



electrode), with respect to which the tip is biased. The counter-electrode is a planar and 
perfect conductor. The tip is also assumed to be a perfect conductor which is biased at VJ,. 
In order to preserve the analyticity of the model, the electrostatic boundary- value problem 
is restricted to a tip with a very simple apex geometry, namely: a hemispherical mesoscopic 
part (radius i? ~ 5 nm) on top of which is half-embedded a small spherical asperity (radius 
Ra ^ R). The contributions of the cantilever and of the macroscopic part of the tip to 
which the apex is connected to are assumed to be negligible. This issue will be justified in 
section [Til A[ The tip-surface distance is denoted z^ (typically a few A). 

The electric field E produced locally between the tip, the dielectric, and the counter- 
electrode polarizes the ionic crystal, which acquires a macroscopic polarization P = n^pi 
oriented along E . In the former equation, n„ is the number of polarizable species per unit 
volume and pi , the local dipolar moment per unit cell. In KPFM, the bias voltage H is 
modulated at frequencies which do not influence the electronic polarization of the ions. The 
major part of the crystal polarization has rather an ionic character, i.e. a net displacement 
in opposite directions of the ions due to their charge ±g with respect to their equilibrium 
positions, ±5/2, with 6 = pi/q, q being the elementary electrical charge. The polarization 
effect occurs as well at the crystal surface, where the positions of the ions become modulated 
perpendicularly to the surface plane, i.e. 5 — > 6-^, as sketched in figsl2]D and d. 6^ is 
proportional to the local electric field Ei at any ionic site and to the total polarizability of 
the dielectric restricted, in our approach, to the ionic polarizability aj (ref.|35|). Usually, Ei 
differs from the external electric field E due to the biased tip because Ei explicitly depends 
on the polarization of the dielectric. The Lorentz's model links Ei and E (ref.[36t]), thus, 
Pi is written: 

^ = aieoEi = XdE = q5-^ — , (2) 

with Xd = «ieo/(l ~ n^ai/S), the dielectric susceptibility of the sample. In the former 
equation, it is important to notice that pi depends on E and hence on the bias voltage VJ,. 
Consequently, this is also true for 6-^. For the sake of discussions, the bias dependence will 
henceforth be explicitly outlined 6-^ — > S-^iVb). Despite S-^{Vb) cannot be estimated at this 
point, it is crucial to keep in mind that the sample surface is polarized by the influence of 
the bias since this is a key aspect of the origin of the CPD atomic contrast. 

Our approach of the electrostatic boundary- value problem relies on an ad hoc assumption. 



The tip, being a perfect conductor, develops a surface charge density a, the origin of which 
is spht into two main contributions a = am + cr^, namely: 

• a charge density am due to the capacitor formed between the biased tip and the 
counter-electrode with the dielectric in between (c/. figlJ]:). Owing to the distance 
between the electrodes, i.e. roughly the dielectric thickness hd, am has a mesoscopic 
character. It is not influenced by the local structure of the tip apex, but rather by its 
overall shape. 

• a charge density a^ originating from the Madelung potential that expands at the crystal 
surface. When the tip is at a distance z^ where the potential is effective, typically a few 
Angstroms, then it develops, in addition to 0"^, a surface charge density a^ {cf. figl2]i). 
a^ has a microscopic character and must strongly depend on the local structure of the 
tip apex and on z^. 

Despite the simple tip geometry that has been assumed, the calculation of the electro- 
static force acting on it due to the combined influence of the capacitive coupling and of the 
Madelung surface potential has no exact analytical solution. Nevertheless, one can build 
up an approximate solution to the boundary-value problem relying on the superposition 
principle. For that purpose, the problem is split up into two boundary-value sub-problems 
schemed in figsJlt and d: 

• problem A: the tip biased at VJ, in front of a dielectric continuous medium (height hd, 
dielectric permittivity e^) held on an infinite planar conductor, the counter-electrode 
(figEb). The local structure of the dielectric is not supposed to influence the tip. 
This is the description of the "capacitive" problem, the solution of which provides the 
surface charge density am- 

• problem B: the tip now biased at V with an infinite plane of alternate point charges 
located at the same distance than the surface of the dielectric in problem A, i.e. z^. 
The layer of point charges is polarized under the electric field that occurs in problem A 
(figl2]l). This describes the "microscopic" problem. The solution provides the surface 
charge density a^. 

Besides, in order to carry out the calculations more easily, it is convenient to distinguish 
the two following geometrical areas on the tip [cf. figl2]i), namely: the asperity, area (1), i.e. 



a half-sphere with radius Ra and the mesoscopic tip apex around it, area (2), a sphere with 
radius R ^ Ra (typically R/ Ra > 50). These two areas do not overlap, but the continuity 
between them is ensured. The vertical force acting on the tip^ is written in terms of am 



and a^ as described above: 



"""' ""^^ -n.uJS = Fm + Fm^ + F^ (3) 



tip 2eo 

n and Uz are the normal to the tip surface and the unitary vector along the vertical z 
axis of the problem, respectively. Doing so, we only focus on the vertical resultant of the 
force acting onto the tip. The above expression can be expanded into three components: a 
purely capacitive part, Fm originating from the tip/dielectric/counter-electrode capacitor; 
a coupling part, Fm^i-, which can be interpreted in terms of the resulting force of all the 
elementary forces due to the electric field (Jm/eo onto each elementary charge a^dS produced 
on the tip by the influence of the Madelung potential of the crystal, Vg] and a purely 
microscopic part, F^, standing for the short-range electrostatic force due to Vg. 

A. Problem A: estimation of (T^ 

Although the boundary-value problem of a planar conductor biased with respect to an- 
other one with an incomplete dielectric layer in between yields an expression of the surface 
charge density, the problem with the sphere does not. However, one can argue that the ex- 
pression of o"m must be a combination between a configuration in which there is no dielectric 
medium in the sphere/counter-electrode interface and an opposite one, where the interface 
is completely filled with it. One can therefore postulate an effective dielectric permittivity 
Q = Kedi where K (< 1) is a constant to be set. Owing to the fact that the mesoscopic 
part of the tip apex, referred to as area (2), is located at a distance Zm = z^ + h^^ R from 
the counter-electrode, the analytical expression of am , explicitly calculated in refs. |37lJ38| . 



asymptotically trends towards the surface charge density of an isolated, biased sphere^^: 

(2) 2m>-R ^d^O^b / -s 

m p V / 

To get the surface charge density on area (1), ain , we seek the potential Vm which equals 
Vb over the asperity. For that purpose, a spherical coordinate system (r, 9, ip) centered on 



the asperity is chosen. The problem having an azimuthal symmetry, the sought potential 
can be expanded in Legendre polynomials^: 

Vl^Hr,^) = f: Unr'' + ^^ ^n(cose) (5) 

n=0 ^ ^ 

At large distance from the asperity, r ^ Ra, the potential must be similar to the one of a 
sphere with radius R biased at Vf,, namely^^I: 

oo 

vJ^\r,e)^>'^>^^v,J2{^yPnicose) (6) 

n=0 

Hence, the coefficients a„ of the expansion are known. The coefficients /5„ are deduced 
from the property of orthogonality of the Legendre polynomials at the boundary condition 
V{r = Ra) = Vfc- The potential is finally written: 



iOO r- 
1+E ©" 



R 



2n+l 



Dn»,n+1 



PnicOsO)} (7) 

n=l '- -• 

The former equation rigorously describes the potential of a system apex/asperity where the 
junction point between the two spheres is smoothed and not singular, as sketched in figl2]i. 
Nevertheless, for R ^ Ra-, the influence of the singular area is negligible. Therefore, equJT] 
is a good approximation to the boundary-value problem. The normal derivation along the 
surface of the asperity yields the expression of o"m : 

~ T r oo / D ^ "^1 



R ^'"^ '' \R 

n=l 



E(2^ + l)^ ^.(cos^) (8) 



Owing to the condition R ^ Ra, the sum can be restricted to the flrst term. Therefore: 

(1) 3ec(eoVf, , , 

"" ^~^° ^^ 

Thus, despite am develops on the asperity, its strength is governed by the radius of the 
mesoscopic part of the tip apex, R and not by the local radius of curvature of the asperity, 

Ra- 



B. Problem B: estimation of cr^ 

The calculation of the surface charge density a^ on the mesoscopic sphere+asperity is 
more tedious, primarily because it relies on the estimation of the Madelung potential of 
the ionic crystal, Vg. The boundary- value problem is now restricted to the determination 
of the surface charge density arising on a metallic tip at zero potential under the influence 
of an infinite planar slab of point charges. Again, the solution of such a problem has no 
straightforward analytical solution. However, we can again use the condition R ^ Ra, as 
depicted, to some extend, in figl2]l. Consequently, the mesoscopic part of the apex can be 
assumed as equivalent to an infinite planar conductor, at least in a small area along the 
sides of the asperity (light grey area in figI2}i). The new boundary- value problem defined by 
an infinite planar conductor influencing another infinite planar conductor at zero potential 
carrying a hemispherical bump with a radius Ra, now yields a solutionis. The method of 
the image charges is used to solve it. The first set of image charges ensuring a zero-potential 
value on the counter-electrode^ produces an electric field which influences the tip. But as 
a matter of fact, this contribution can be neglected because the distance between the tip 
and the counter-electrode is on the millimeter range and the image charges originate from 
the Madelung potential, which is known to decay exponentially fast^i^ (c/. also hereafter). 
The second set of image charges^ is quasi-punctual and located at the center of the asperity. 
Hence, owing to the simplified geometry of the electrode, the problem is reduced to a sphere 
with radius Ra at zero potential in the influence of two infinite planes, i.e. the slab and 
its image, which are anti-symmetrically spaced with respect to the sphere. This procedure 
ensures a zero-potential on the approximated plane within which the asperity is embedded. 
Again, it is more convenient to use a spherical coordinate system centered on the asperity. 
Then, the surface charge density a^ = al +al is derived by normal derivation along areas 
(1) and (2), namely: 



cr„ 



eo- 



dV,{r,e,v) 
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(l)=aspcrity 
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dV^{r,d,^) 
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r=Ra 



{2)=planar area 



(10) 



0=7r/2 



The potential V^ is derived from the Madelung potential of the ionic crystal, previously 
referred to as K, by the method of the image charges and the Kelvin transform (influence 
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on a sphere biased at zero potential)^, namely: 



V,{r, 9, v) = [Vs{r, 9, ^) - ^V, (^, 9, <^) } 
|K(r, 7r-9,^)-^Vs(^^,7r- 9, y.) } 



slab 



(11) 



image slab 

The former equation fulfills the boundary condition V^(r = Ra) = everywhere along the 
surface asperity or along the surface of the local planar area around it. 

Vs can be estimated on the base of the work by Watson et al.^^. When considering an 
infinite planar slab of point charges, the authors state that the potential, so-called Madelung 
surface potential, reaches its asymptotic value in a very short distance normal to the slab. 
Consequently, the ions within the crystal at a distance only one lattice constant from the 
surface have Madelung potentials which are indistinguishable from those of the bulk. In 
other words, the tip will mainly be influenced by the surface potential and not by the one 
arising from the bulk part of the ionic crystal. This is why a single, infinite, layer of point 
charges is enough to describe the infiuence of the Madelung surface potential on the tip, 
which motivates our initial assumption. The potential is written: 

where p = xi + yj is the polar vector of any ion of the surface slab in an orthogonal basis 
{0,i,j), O being the projection of the center of the asperity on the surface, i and j the 
unitary vectors of the fee unit cell. The summation is performed over the vectors G of 
the reciprocal lattice of an arbitrarily defined unit cell and a' is proportional to the lattice 
constant a of the fee unit cell. q{G) is a structure factor: 

qiG) = ^J2^,e^-'^^^^'''^ (13) 

k 

It is summed over the ions within the defined unit cell. The k*^ ion carries an electrical 
charge g^. Its planar and perpendicular coordinates from the origin of the cell are given by 
the two vectors S I and 5 ^. The latter refiects the polarization effect felt by the ion within 
the unit cell, previously referred to as 6^{Vb). In equITSl it is assumed that the polarization 
of the ions extends all over the (001) surface plane. In any case, it must extend over an 
much larger area than the tip asperity. This assumption is consistent with the electric field 
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produced by area (2), which is constant within an area roughly scahng as the niesoscopic 
tip radius R. Vg is calculated from the unit cell defined in fig|2^ (light grey). It consists of 
4 anions and a cation weighting for a fourth and one, respectively. The vectors of the direct 
lattice are 7^ = a'i and /3 = a'j, where a' = a^/2/2. Owing to the exponential decay of the 
potential with z^, visible in equJT2l the calculation of the structure factor can be restricted 
to the first four reciprocal vectors, namely: Gf^^. ^ = ±27r/a'(i or j). The calculation yields: 

Vsix, y, z^) = ^ cosh[3'^(H)]x(a;, y)e-^'^ (14) 

with: 6^(yb) = ^S^iVh) and xi^jV) = cos [^{x — Xq)] + cos [^{y — yo)], a spatial modu- 
lation term, xq and yo are the x and y coordinates of the center of the asperity projected 
onto the unit cell. Setting xq = yo = locates the asperity and therefore the tip on top of 
an anion, the reference ion within the defined unit cell. The above expression exhibits the 
expected exponential decaying behavior as a function of z^. The potential is reported in 
figj3]for a = 0.66 nm, 6-^ = 11 pm and z^ = 4 A. The value of 6^ will be justified in section 
lYl 

The expression of a^ = cr}t + 0")^ can now be derived from equilOl The calculation of 
all yields: 

C(i?.,^,^)^(^)(e)}e-!^(^-+^") 

with x(r, 0,^9), the expression of x(x, y) in the spherical coordinate system centered on the 
asperity and ({r, 6, ip), the spherical expression of another spatial modulation term given by: 



C{x,y) = (xsin [^{x — Xq)] +ysm [^{y — yo)])/\/x'^ + y"^. The functions JF^^) and ^"^^^ are 
written: 

2a' . , fr]e\ a n ^ fVe 



jrm (^) = ^ ginh (11) - 8 COS ^ cosh i"^) (16) 

vrila V 2 / V 2 / 



?(l)('fl^ = SclnflclnV, (^ 



and: 



^^^)(^) = 8sin^sinh(^j (17) 

with: Tje = Arc Ra cos 6 /a' (see also the appendix). Thus, it was necessary to assume area (2) 
as an infinite plane in order to derive the expression of aji by means of the method of the 
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images. In such a case, the expression of cr}, is rigorously derived from equ{TO] and is exphc- 



itly given in ref.[43(]. For x and y positions large compared to Ra however, this description 
does not fit with the geometry of the tip apex defined in problem B. The "infinite plane" 
must actually be shrunk down to a spatially limited area around the asperity, as sketched 
in figEli. This is made possible when assuming that any planar area is the asymptotic limit 
of a sphere with large radius compared to its extension. The former statement is fulfilled 
by the condition R ^ Ra- The surface charge density of a sphere with radius R under the 
influence of the surface potential Vs is also derivable from the Kelvin transform. Here, the 

(2) 

spherical coordinate system is centered on the sphere with radius R. a}j, is written: 

af = ^cosh[5^(H)] {r^'KomR^e,^)- 



The functions T^"^' and Q^"^' are written: 



a . „\ 27r : 



TxR 



and: 



^(2)(^) = 4sin0e-^«(™^^+i) (20) 

Since !F^'^\0) and Q^'^\0) decrease exponentially fast as one moves away from the foremost 
position of the tip apex, it can readily be verified that equHH] and the expression given in 
ref.[43| do fit for x oi y E \R„:,'iRn\. Thus, even in the vicinity of the asperity, expression [T8l 



can be used instead of ref . 43|] . Therefore, the most part of the contribution of the sphere 
to ail is restricted to a small area that can be assumed as locally planar. 

The graph of the projection of o"^ on the tip apex (areas (1) and (2)) is reported in figH] 
on top of an anion (positive charge density on the asperity) at a distance z^ = A A. from the 
surface and for R = 5 nm and -R^ = 1 A. The oscillations of afj_ at the surface of the tip due 
to the image charges of the crystal are readily visible, although their amplitude decreases 
exponentially fast along the sides of the asperity. The relevant oscillations in area (2) spread 
out in a circular area with only two unit cells radius. Owing to the exponential decay of the 
Madelung surface potential farther away from the asperity, the tip surface is too far from 
the crystal surface to produce relevant image charges. On top of the asperity, the surface 
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charge density reaches about 2 x 10 ^ Cm ^. For a half-sphere with radius Ra = 1 A and 
with (5-*" = 11 pni, this results in a local charge carried by the tip of about 0.8 x 10^^ q. 



III. ESTIMATION OF THE FORCE 



A. Estimation of Fm 



Owing to the geometry of the problem and since am is equivalent to the surface charge 
density of an isolated, conducting sphere, we can now fully evaluate the purely capacitive 
component F^ to the total force. The force effectively acting on the mesoscopic part of the 
tip apex can be derived from the image charge of the sphere —AirR'^am placed at a symmetric 
position with respect to the counter-electrode, i.e. at a distance 2hd from the mesoscopic 
apex: 



(47ri?V„)2 ttR 



2 



47reo(2/irf)2 h 



--^.eoV,' (21) 



With ea ~ 4.87 (ref.|44]) and K = 0.9 (arbitrarily), then e^ = Kea = 4.38. For a typical 
ac modulation of the bias of about 0.5 V and considering i? ~ 5 nm and h^ = 5 mm, the 
former equation gives an estimation for the mesoscopic component of the electrostatic force: 
Fm ~ 8 X 10~^^ pN. As suspected, this contribution is negligible compared to the other two 
components^. Therefore the expression of the total force simplifies to: 

F = Fm, + F,= [ ^^^dS+ [ ^dS (22) 

J tip ^0 J tip ^*^o 

The reason why the tip has been restricted to its apex is now clear. Regarding the for- 
mer capacitive force, owing to the distance between the tip and the counter-electrode, the 
contribution of a macroscopic body in addition to the mesoscopic apex would not change 
notably the total force. Regarding the other force components, Fm^ and F^, owing to the 
exponential decay of the Madelung surface potential, the influence of the macroscopic body 
of the tip is expected to be negligible as well. 
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B. Estimation of F^ 

Let us first focus on F^. Tlie integral must be performed over the asperity and the 
mesoscopic sphere around it, areas (1) and (2), respectively. It is recalled that the two areas 
do not overlap each together. For each of them, the spherical coordinate system must be 
centered on the corresponding sphere. Thus: 



F^ = FJ') + F^'^ = 1^ /: /q'" cos e sin ^^'^ d9d^+ 
g /!"'' C cos sin Oaj^^^dOdip 
The radial coordinates are fixed to r = Ra and r = i? in areas (1) and (2), respectively. 
The polar integration in area (1) is performed with 6 G [7i/2;7i], whereas in area (2), it is 
performed with 6 G [vr/2; 6m], where 6m = tt — arcsin(i?a/i?), which ensures the continuity 
from (1) to (2). Regarding area (2), the choice of the beginning angle of the interval (7r/2) 
does not infiuence notably the result of the integration. In other words, owing to the 

(2) 

exponential decay of a/j , the exact shape of the mesoscopic part of the tip apex far from 
the asperity is not relevant. This justifies a posteriori the choice, although simple, of a 
spherical geometry. Both integrations over the azimuthal angle must be performed over 
V? G [0; 2tt]. Note that in each of the former integrals, the term cos^ stands for the vertical 
projection of the force, as stated initially. Although the integration over the azimuthal angle 
yield an analytical result, the integration over 6 does not, which requires to evaluate some 
integrals numerically. The expression of F^ is written: 

^i'^ = fiS cosh^[^^(H)]e-^(^''+«") {A(i)+ ^^^^ 

B(^) [cos (xo) + COS (yo)] + 2C^^'> cos (xq) cos (yo)} 
with: Xq = 211x0/ a' and y^ = 27ryo/a', the reduced coordinates of the tip above the crystal 
surface. The integral forms of coefficients A^^\ B^^' and C^^' are reported in the appendix 
as functions of Ra and a' . Taking a typical lattice constant for alkali halides a' = a^/2/2 ~ 
0.45 nm and assuming Ra — 1 A, we get: 

A(^) ~ -130 5(^) ^ -70 C(^) ~ A(^) ~ -130 (25) 

Thus, F^ explicitly depends on the spatial modulation of the surface potential. Note also 
the exponential decay with the distance, actually faster than the distance dependence of the 
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Madelung surface potential and also the doubling spatial period term. 
Similar integration on area (2) yields: 

Fl^ = fiS cosh^[?^(H)]e-^(^''^«") {A(2)+ ^^^^ 

5(^) [cos (xo) + COS (yo)] + 20'^^^ cos (xq) cos (yo)} 

The integral forms of coefficients A^'^^ , B'^'^^ and C-^^ are derived from those of coefficients A^^^ , 
fi(^) and C(i), by replacing Ra, J^^^\e) and g'^^\e) with i?, J^(2)(^) ^nd G^'^^e), respectively. 
The integration is now performed with 9 G [7r/2; 9m]- With similar parameters than before 
and setting R = 5 nm, we now get: 

A(2) ^ _8 5(2) ~ C(') ~ (27) 

The striking discrepancy between A^'^\ B^"^^ and C*-^-' is due to the combined contribution of 
the Bessel functions and the exponential decay of the functions J-'^'^\9) and G^'^\9) occurring 
in the integral forms of the coefficients B^"^^ and C^'^\ Hence, the force simplifies to: 



ir^2) ^ Q_^ cosh2[5^(H)]e-!^(^''+^")A(2) (28) 



2eoa'- 

The spatial modulation of the potential does not influence the mesoscopic part of tip while 
scanning the surface. This contribution acts as a static shift to the total electrostatic force, 
similarly as the Van der Waals long-range interaction for the short-range chemical interac- 
tions which are responsible for the topographic atomic contrast in nc-AFM. 

C. Estimation of Fm^ 

The geometrical splitting in terms of areas (1) and (2) used for the estimation of F^, can 
equivalently be applied to Fm^. Thus: 

p _ p(l) I p(2) _ 

^ j: J^^ COS 9 sin 9a^^^ai}''d9d^+ (29) 

f I^"" C cos 9 sin 9ain^ai^^d9d<f 

am and am are the surface charge densities on areas (1) and (2) within the frame of problem 
A. The calculation of Fm/i yields: 
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L)(^) [cos (xq) + cos (yo)] 

The integral form of D^^' is reported in the appendix. The numerical integration results in 
D^^^ ^ —15. Similar integration over area (2) gives a similar expression of Fm/i by replacing 
indexes (1) with indexes (2) and Ra with R. The integral form of the coefficient D^"^^ is 
similar to the one of D^^\ except that the term cos^ 6' is replaced by cos 6'. The numerical 
integration of -D*^^-* yields almost zero. This was expected since the integration is performed 

(2) 

over the oscillations of the charge density a^ due to the image charges of the surface 
potential (c/. fig SI). On the opposite, this was not observed in F/i because the integration 
was performed on the square oi ali . Therefore Fm^ is finally written: 

F = F^^^ (31) 

Hence, the former equation states that the coupling of the force to the bias H actually 
occurs only by means of the microscopic surface charge density at the foremost part of the 
tip, i.e. on the asperity. Owing to the integration over the mesoscopic part of the tip apex 

(2) 

and the subsequent cancellation due to the oscillations of a/t , no relevant coupling between 
the bias voltage and the mesoscopic part of the tip apex can occur. This aspect strongly 
suggests that the effect we are reporting is mainly controlled by the foremost structure of 
the tip. The effect is expected to be much enhanced in case of tips with sharper geometries, 
particularly those with apexes including atomically sharp edges. 

The expression for the vertical contribution of the total force acting on the tip due to 
the combined influence of the capacitive coupling and of the Madelung surface potential is 
flnally written: F = F^^ + F^ = F^ + F^^^ + F^^^ with F^ given by equISOl F^^^ by equl2l 
and F^ by equj28l The graph of the total force and of its components is reported in flgl5] 
for similar parameters than previously, namely: a = 0.66 nm, 5-^ = 11 pm and z^ = A A. 
A typical value of bias has been set, Vb = +1 V. Thus, it is visible that the term F^ is 
negligible compared to others. This is due to the prefactor i?^/a'^ vs. R^/a"^ for F^ . The 
total electrostatic force F flnally simplifles to: 

F = F^l + Ff^ (32) 
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In figJSl the force reaches an average value of about 9 pN (absolute value) and a corrugation 
of about 2 pN (peak to peak). 

IV. IMPLICATIONS FOR KPFM 

A. Estimation of 5^ (Vfc) 

The net displacement S'^iVb) of the topmost ionic layer induced by the polarization can 
be estimated out of the electric field E between the tip and the surface (c/. equH]). Section 
III Al has shown that the field induced by the bias voltage was essentially controlled by the 
mesoscopic radius of the tip, R, and not by the foremost asperity. Therefore, owing to equsH] 
or[8j -E ~ Vb/R. The expression of S^iVb) can now be deduced using equj2j 

6^{Vb) = ^E = VM (33) 

q qR 

An order of magnitude for S-^{Vb) can now be calculated as follows: with Vb ^ IV, R = 5 nm, 

a, = 70 X 10^'^° m^ (ref. 36|) and n^ = 8^/2/ a^ ~ 40 x 10^^ m~^ (number of polarizable ionic 

species per volume unit in a fee crystal of KBr with a lattice constant a = 0.66 nm), we get 

E = 2 X 10^ V.m"^, Xrf ^ 9 X 10"^^ F.m^ and therefore 6-^ ~ 11 pm. 

B. Detected signal in KPFM: connection with the local CPD 

When performing KPFM experiments, the bias voltage Vb is modulated with a frequency 
/fc and may as well include a static component to compensate for the long-range electrostatic 
forces. This is the reason why the average value of the experimental CPD image shown in 
figil^ reaches -3.9 V. It is not rare that, on ionic surfaces, many volts are required to 
compensate for the long-range electrostatic forces due to trapped charges while the cleavage 
of the cristal^^. Thus: 

Vb = Vd, + VacSm{27Tfkt) (34) 

The electrostatic force is thus triggered at fk and then detected as an additional low- or high- 
frequency component when doing FM- or AM-KPFM, respectively. In both techniques, a 

(c) 

proper dc bias voltage produced by an external controller, hereafter referred to as Vj^ , is 
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applied between the tip and the counter-electrode to cancel the modulated component at 
fk, i.e. the oscillation amplitude of the second bending eigenmode of the cantilever in AM- 
KPFM, or the one of the frequency shift in FM-KPFM. When applied to the tip {i.e. with 
the counter-electrode grounded), this dc bias is the opposite of the local CPD defined in 
equHJ Vj^'' = -Kpd. 

In order to stick to the AM-KPFM experiments (fig{T]), it is necessary to estimate the 
amplitude of the second eigenmode of the cantilever and then derive a condition on the dc 
value of the bias ultimately able to nullify it. On the one hand, it is therefore mandatory 
to check carefully all the occurrences of the modulated component of the bias voltage in 
the expression of the force. This includes explicit dependencies, such as those due to the 
polarization, but also implicit ones, as discussed hereafter. On the other hand, an expression 
of the oscillation amplitude of the mode modulated at fk must be derived. 

We first address the problem of explicit and implicit bias dependencies in the expression 
of the electrostatic force. Owing to the explicit VJ, dependence in 6-^, which has been kept 
throughout the description of the model, it can be seen that the polarization effect is mainly 
included in the cosh function of Fm/i and F^ through a linear and a quadratic dependence, 
respectively (c/. equs J30] and [281) . Furthermore, since 6'^ is small compared to a', cosh[(5^(Vfe)] 
can be expanded in series. To first order: 

cosh[S\V,)] ^ 1 + 5\V,)' = 1 + ^S^iV,)' = 1 + ixM' (35) 

with x'd = '^'^Xd/{0''<lR)- Replacing this expansion in the expressions of the components of 
the force and keeping the linear and quadratic terms in Vb yields, with compact notations: 

i^s=?'.i^s*a^ (36) 

and: 

i^f =^f*?^[l + 2(x'.H)^]^, (37) 



where Kmli and Ki^ are two dimensionless coefficients standing for geometrical factors of 
areas (1) and (2), respectively: 

A-S = ^C'" aud /<-f = ^A« (38) 
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$m/i, $}j are also two dimensionless coefficients carrying the spatial dependence of each 
force component: 



$S = e-^(^-+^") [cos (xo) + cos (yo)] (39) 

and: 



$(f) = e-^(^M+«.) (40) 

Implicit 14 dependencies are now discussed. In the two former equations, particular attention 
must be paid to z^. So far, this parameter was defined as the tip-surface distance and set to 
an arbitrary, constant, value. However, when dealing with AM-KPFM, 2;^ is not static but 
actually coupled to the bias and to the oscillation amplitude of the fundamental eigenmode 
of the cantilever. In the following, for the sake of clarity, the variables related to the 
fundamental eigenmode of the cantilever will be denoted with index "0" and those of the 
second eigenmode with index "1" . Thus, let Zq, Z\ and D be the instantaneous position of the 
fundamental eigenmode of the cantilever, the instantaneous position of the second eigenmode 
of the cantilever and the distance between the surface and the equilibrium position of the 
cantilever at rest, respectively. Therefore, z^{f) = D — zo{t) — zi(t). Hence, if Vb has the 
form given in equlMl one can postulate, to ffist order: Zi{t) = Ax sin(27r/fct + ty^i). A^ and 
V?i stand for the oscillation amplitude of this mode and its phase lag with respect to the 
electrostatic actuation, respectively. Their exact expressions are not easily derivable, but it 
must be noticed that A\ must be connected to the amplitude of the modulation, namely 
Vac- When jy, accurately matches the actual resonance frequency of the second eigenmode^, 
then (y9i = —-njl. zq is experimentally driven by the control electronics of the microscope at 
the actual resonance frequency of the fundamental mode of the cantilever, /q. It is known 
that it has an almost harmonic behavior of the form: Zo{t) = Aq sin(27r/ot — vr/2). Thus, 
$m/i and $/! have a component at /o which is further modulated by the dynamics of the 
second eigenmode, electrostatically actuated at fk- 

We can now propose a self-consistent approximated solution to the equation of motion 
for Zi{t) and thus derive the expression of the oscillation amplitude Ai. This equation has 
the standard form: 
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^■i(t) + ^^lit) + ^'Mt) = ^ + £ii!t±Ij^^ (41) 

Qi mi nil 

where Fext is an external force oscillating at /o which controls the dynamics of zo{t). Let us 
assume: i- Ai <^ Aq, i.e. the amplitude of the mode is much smaller than the one of the 
fundamental mode, ii- Ai <C a', where a' = ay2/2, a being the lattice constant of the crystal 
and iii- that the dynamics of zi{t) is mainly influenced by components at fk- Assumptions 
i- and ii- are not too strong, since the experimental estimations of Ai yield a few tens of 
picometers. Assumption i- implies that the dynamics of the fundamental mode is not much 
influenced by the one of the second eigenmode. Hence, the solution of the equation of 
motion of ZQ{t) has indeed the form postulated above. Assumption ii- allows us to linearize 
$m/i and $/i with respect to zi(t). Finally, assumption iii-, which is consistent with the 
postulated solution for Zi(t), zi(t) = Ai sin(27r/fct+y9i), simplifies further the above equation 
of motion. Now, owing to the former assumptions, equJH] can be solved by injecting the 
postulated expressions of both eigenmodes and keeping only the terms oscillating at fk- For 
that purpose, the exponential term wherein zo{t) occurs must be expanded in Fourier series. 
Then, the only possibility to preserve terms at fk is to keep the lone static component of the 
Fourier expansion, hereafter referred to as oq (expansion of ^mfj,) and 60 (exansion of $}j ). 
After linearization, $m/i and $}j can finally be written as: 



1-, 11- 

00 



{ao + ^ an cos{2TTnfot) } [cos (xq) + cos (yo)] ^ ^ 



n=l 



neglected, owing to iii- 



and: 



$(2) ^ 



— ^iW 



g-^(D-Ao+H„)|^^ ^ ^ ^^ cos(27™/ot)} (43) 



n=\ 



neglected, owing to iii- 



With: 
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/ is the modified Bessel function of the first kind. When replacing equsJl2] and |13] with 
zi{t) = y4i sin(27r/fct + ipi) in equE], it is possible to derive an expression for Ai. The 
condition on V^^ to match Ai = is finally written: 

Vr = -^f^^lie-"'^-'"'' |cos(f„) + cos(,-„)| (45) 

The graph of V^^ is reported in fig|6] for a = 0.66 nm, S^ = II pm, Aq = 5 nm (hence 
oq = 0.0487 and Bq = 0.0344) and {D — Aq) = 3.5 A. The value of Aq has been chosen 
consistently with the experimental conditions. On the contrary, the tip-surface distance 
has been chosen arbitrary, but however in a range where the atomic contrast is usually 
experimentally achieved. In figEl the lateral periodicity of the underlying lattice is readily 
visible, but surprisingly, the potential scales between -0.6 to +0.6 V from an anionic to a 
cationic site, respectively. At similar height, this is three times larger than the Madelung 
surface potential (c/. figint). The comparison with the experimental results is more severe 
since the theoretical prediction is one order of magnitude larger. At this point, it is recalled 
that the strong tip geometry dependence of the problem makes a straightforward comparison 
between the theoretical prediction (equ HSl) and the experimental results difficult, since our 
analytical expression of VJ^ = — Kpd relies on a somewhat unrealistic tip. 

C. Experimental implications for AM-KPFM 

The figure and the above formula show that the bias voltage to be applied on the tip to 
compensate for the electrostatic force is governed by three main factors: 

• the dielectric properties of the sample such as its dielectric permittivity and lattice 
constant. 

• a subtle balance between mesoscopic and microscopic geometric factors of the tip. 

• a lateral periodicity similar to the Madelung surface potential of the crystal. 

A straightforward consequence is that the atomic corrugation of the CPD reported ex- 
perimentally might stand for the spatial fiuctuations of the Madelung surface potential, 
however with an amplitude that depends on the surface polarization and hence on the ap- 
plied ac voltage. It is also important to notice that the CPD compensation is proportional 

22 



to Kml, i.e. to the asperity size. The former being also the source of the couphng between 
the tip/dieletric/counter-electrode capacitor and the Madelung surface potential, i.e. the 
source of the KPFM signal, the atomic contrast of the CPD is therefore closely connected 
with the geometry of the very foremost part of the tip. This is consistent with the short- 
range character of the interaction. But, on the other hand, the explicit dependence with 
geometric factors of the tip, unambiguously proves that quantitative measurements of the 
local CPD are unlikely to be performed in KPFM, unless the tip geometry be accurately 
known, which is practically never true. 

Furthermore, although equJ45] explicitly exhibit a distance dependence, consistently with 
the experimental observations in the above mentioned references, an increase of the compen- 
sated CPD as a function of the distance is nevertheless surprising. The residual exponential 

(2) 

dependence originates from $/( , i.e. from the influence of the Madelung surface potential on 
the mesoscopic tip apex. As discussed above, the CPD compensation being partly governed 
by the asperity, a tip-surface distance increase Az^ produces a decrease of the related force 
Fm/i proportional to exp {—2TTAz^/a') (of. $m^, equsJ39]or |42|) . This abrupt change is com- 
pensated by an equivalent exponential increase of the compensated CPD. This process can 
obviously not occur at any tip-surface distance. Prior to being cancelled, the fk component 
must stand for a measurable signal. Therefore the above discussion stands within a narrow 
range of distances from the surface, typically a few times the asperity radius. 

Nevertheless, an increase of the measured local CPD as a function of the distance is still 
expected to occur if the distance dependence in $}i decays faster than the one in ^mii- In 

(c) 

other words V^^ must increase with the distance as soon as the distance dependence of the 
force induced by the influence of the Madelung surface potential on the mesoscopic part of 
the tip decays faster than the one of the force induced on the asperity due to the capacitive 
coupling with the image charges of the surface. 

Let us finally point out that such a distance dependence of the CPD might make the 
experimental achievement of the atomic contrast easier, which strengthens the argument of 
an intrinsic imaging process of the local CPD. Indeed, no major tip and/or surface distortion 



is expected to occur in an equivalent range of distance (> 3 A, ref. J47l]). In that case, 
instabilities due to adsorbed and/or mobile atomic or ionic species at the tip apex are less 
likely to occur, which makes the imaging process robust, as experimentally observed. 
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To conclude, the analytical approach, although restricted to a tip with a basic geometry, 
remains helpful, primarily because it provides an expression of the short-range electrostatic 
force that can be connected to the nc-AFM-KPFM simulator. There are obvious limitations 
to our approach, the most important one being the use of classical, continuous electrostatics 
to treat the angstrom-size nanoasperity. This obviously must break down at a certain point, 
and be replaced by a proper quantum mechanical treatment of the problem. In the near 
future, the electrostatic model should be extended to a bit more complex systems such as 
local dipoles, charges or defects at the surface and at steps of ionic crystals. 

V. CONCLUSION 

The aim of this work was to provide a consistent approach to describe the short-range 
electrostatic force between the tip of an nc-AFM and the (001) surface of a perfect ionic 
crystal. In order to develop an analytical expression for the total electrostatic force, the tip 
has been restricted to a simple geometry and the influence of the sample has been described 
by means of its Madelung surface potential. In such a way, an analytical solution for the total 
electrostatic force was found within the boundary- value problem assuming a thick dielectric 
sample and an infinite top-layer of ionic surface charges. 

Two major contributions to the electrostatic force can be extracted: the first stands 
for a coupling term between the microscopic structure of the tip apex and the capacitor 
formed between the tip, the dielectric ionic crystal and the counter-electrode due to the 
bias voltage VJ,; the second term depicts the influence of the fluctuations of the Madelung 
surface potential arising at the surface of the ionic crystal on the mesoscopic part of the tip, 
independently from its microscopic structure. The former has the lateral periodicity of the 
Madelung surface potential whereas the latter only acts as a static component which shifts 
the total force. 

Beyond the dielectric properties of the crystal, which are explicitly included in the model, 
the ionic polarization of the sample due to the influence of the bias voltage applied to the 
tip/counter-electrode capacitor is mainly responsible for the atomic contrast of the KPFM 
signal. Typical orders of magnitude give a net displacement of the ions of about ±10 pm 
from their equilibrium positions. Note that this displacement only occurs if a tip-sample 
bias is applied (ac or dc), which is always the case in KPFM experiments. 
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A detailed analysis of the bias voltage required to compensate for the electrostatic force 
shows that the compensated CPD has the lateral periodicity of the Madelung surface po- 
tential. However, there is a strong dependence on the tip geometry, the applied modulation 
voltage as well as the tip-sample distance, which can even lead to an overestimation of the 
real surface potential. 

For a quantitative evaluation of KPFM results, it is thus essential to account for all 
the parameters of the experiment, among which the tip shape. The analytical expression 
developed in this work provides an alternative tool to elucidate the contrast formation in 
KPFM on ionic crystals, and in combination with the nc-AFM simulator it might enable us 
to interpret our results more accurately. 

APPENDIX 

After integration of F^ over the azimuthal angle cp E [0; 27r] , we have: 

F« = g|cosh2[J±(y^)]e-l^(^.+i?.) 1^(1) + 
5(^) [cos (xq) + cos (jjo)] + 2C(^) cos (xq) cos (yo)} 



(46) 



where: 



/'TT 

A^^^ = TT / COS ^ sin ^[2J^(^)' (6) + G^'^^' {e)]de (47) 

^1 



5(1) = n J^ cos ^ sin ^ { (^(^^ (9) - ^^) Jo (Ve) - 
2^«(^)^«(^) Ji (r^e) + ^^J2 ive)] dO 



(48) 



and: 



C« = 7r/:cos^sin^{2J^W'(^)Jo(r7^)- 

2v^^«(^)^«(^)Ji(^^)+ (49) 

with: ri0 = AirRaSinO/a' and r]g = 2\/2-KRaSm9/a'. Jo, Ji and J2 are the Bessel functions 
of the first kind. 

Regarding Fm/i, the integration over the azimuthal angle (^ G [0; 2%] yields: 
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"" - '|l#<:osh[«x(Vi,)|ViB<i) [cos (f rc„) + 



m/i 

COS 



(50) 



where: 



D(i) = -27r/:cos2^sin^{j^«(^)Jo (f ) - 



(51) 
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^^ The latter is obtained by anti-symmetrical transform of the charges layer with respect to the 
counter-electrode symmetry plane. 

^^ R. Watson, J. Davenport, M. Perlam, and T. Sham., Phys. Rev. B 24, 1791 (1981). 

'^^ This is also true because it has been assumed that the ionic crystal carried no net charge. 
Therefore there is no long-range electrostatic interaction. 

^■^ The latter is now obtained by anti-symmetrical transform of the layer of ionic charges into the 
asperity. 

^^ The expression oi a\i {r > Ra,0 = t:/2,lp) for an infinite planar geometry is: ajj, = 
^ cosh[5^(T4)]{4x(i?, vr/2, ip) - ^^xiR^r, vr/2, (/,)}. 

'^^ C. Andeen, J. Fontanella, and D. Schuele, Phys. Rev. B 2, 5068 (1970). 

^^ As a matter of fact, h^ being very large, it can readily be verified that this result is always true 
whatever is the exact shape of the mesoscopic part of the tip apex. 

^^ It must be noticed that the resonance frequency of that mode might be shifted owing to the in- 
fluence of the overall attractive force, i.e. including Van der Waals, short-range and electrostatic 
components, similarly to the resonance frequency of the fundamental mode. 

^"^ In ref.[21], the foremost unstable atom of the tip which is responsible for the contrast in dis- 
sipation is triggerred at a distance of about 3 A from the surface. Above that threshold, no 
instability is observed. 
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FIG. 1: a- Experimental image showing the atomic contrast of the compensated CPD measured 

on a (001) surface of KBr in ultra-high vacuum by means of AM-KPFM. The vertical contrast 
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ranges from -3.95 to -3.85 V from dark to white spots. The dashed line depicts the cross section 
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FIG. 2: a- KBr lattice with a fee structure eorresponding to the dotted square shown in figH^. The 
white spots of the experimental image have been placed on top of anions arbitrarily, b- Scheme 
of the KPFM experiment defining the electrostatic boundary-value problem to be solved. The 
metallic tip is biased with respect to a metallic counter-electrode placed a few millimeters far from 
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FIG. 3: (Color online), a- Madelung surface potential calculated from equll4l for a = 0.66 nm, 
J-*- = 11 pm and z^ = 4 A. The vertical contrast ranges from —100 (blued spots) to +100 mV 
(reddish spots). The unit cell depicted with a dotted line is centered on a cation, b- Distance 
dependence of the potential on top of an anion (dotted curve) and on top of a cation (continuous 
curve) showing the exponential decay of the potential, c- Cross section along the dotted line shown 
in a-. 
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FIG. 4: (Color online) a- Perspective view of the projection of the microscopic surface charge 
density a^ on the tip calculated on top of an anion from equs{T5] and [18] with Ra = 1 A and 
R = 5 nm. a^ is strongly increased on top of the asperity. Then, owing to the exponential decay 
of the Madelung surface potential, the surface charge density strongly decreases, which makes the 
contribution of the mesoscopic apex weak, b- Top view of o"^. The oscillations due to the image 
charges of the ions at the surface of the crystal are well visible. The most part of the attenuation 
of o"^ occurs within a single unit cell. Two unit cells apart from the asperity, a^ is almost zero. c-. 
Cross section along the dotted white line shown in b. 



32 




03 
O 



CD 

-I— ' 

.O 





-1 

-2 
-3 




2, (nm) 



0.3 0.4 0.5 0.6 0.7 0.8 



cross section (nm) 



CD 

o 



.0 




-2 
-4 

-8 



.„f^... 



^V,; 



v^; ;;; 



.^'^ 



-0.6 -0.4 -0.2 0.2 0.4 0.6 

cross section (nm) 



FIG. 5: (Color online), a- Total electrostatic force over a unit cell calculated from F = Fm^ + 
Ffi ' -\- FJi ' for a = 0.66 nm, 6-^ = 11 pm, z^ = 4 A and Vb = +1 V. The unit cell depicted with a 
dotted line is centered on a cation. The vertical contrast ranges from -10 (blued spots) to -8 pN 
(reddish spots). The force is more repulsive on top of cations (central ion) than on top of anions, 
consistently with the bias polarity, b- z^ dependence of the electrostatic force on top of an anion 
(dotted curve) and on top of a cation (continuous curve), c- Cross section along the dotted line 
shown in a- showing the total force (thick continuous line), Fm/i (dashed line), F^ (dotted line) 
and Ff^ (greyed line). F^ is negligible compared to the others. 
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FIG. 6: (Color online), a- V^^ bias voltage required to compensate the local CPD calculated from 
equl45lfor a = 0.66 nm, 5 = 11 pm, Aq = 5 nm and z^ = 3.5 A. The vertical contrast ranges from 
-0.6 (blued spots) to +0.6 V (reddish spots). The unit cell depicted with a dotted line is centered 
on a cation, b- Distance dependence of the potential on top of an anion (dashed curve) and on top 
of a cation (continuous curve), c- Cross section along the dotted line in a-. 
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